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Abstract 



Along the recently trodden path of studying certain number theoretic prop- 
erties of gauge theories, especially supersymmetric theories whose vacuum mani- 
folds are non-trivial, we investigate Ihara's Graph Zeta Function for large classes 
of quiver theories and periodic tilings by bi-partite graphs. In particular, we ex- 
amine issues such as the spectra of the adjacency and whether the gauge theory 
satisfies the strong and weak versions of the graph theoretical analogue of the 
Riemann Hypothesis. 
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1 Introduction 



Though the geometry of gauge theories has constituted one of the greatest intellectual 
triumphs and served as a heart- warming matrimony between physics and mathematics, 
the infusion of number theory, that Queen of mathematics, into or her emergence 
from physics has been comparatively humbl^ Yet, there has been a colourful skein 
threading the two seemingly disparate enterprises since at least Hilbert and Polya in 
addressing the Queen of number theory, the Riemann Hypothesis: how one might find 
a Hermitian operator whose eigenvalues coincide with the critical zeros of the Riemann 
zeta function. 

Over the decades, there have been various proposals in finding analogues of the 
zeta function which satisfy a functional equation symmetric about the vertical line 
situated at ^ in the complex plane, as well as possess suggestive zero free regions. 
The most remarkable of these is indubitably the Hasse-Weil zeta function which is an 
(exponentiated) generating function for algebraic varieties over finite number fields. 
The relation of its functional form to the geometrical data has been one of the greatest 
feats in modern mathematics. Recently, there has been some activity in investigating 
the physics of this zeta function, for compact Calabi-Yau manifolds [l], as well as for 
non-compact vacuum varieties in relation to gauge theory [2]. 

This latter case of supersymmetric gauge theories with non-trivial vacuum promises 
a particularly enticing ground for the inter-play amongst physics, algebraic geometry, 
combinatorics and number theory. Physically, they are generically theories with prod- 
uct gauge groups and comprise of such important examples as the Standard Model 
and those arising from the AdS/CFT Correspondence. A so-called plethystic pro- 



gramme 18 , exploiting algebro-geometric and combinatorial techniques, has been de- 
veloped to enumerate the gauge invariants in such theories, as too was the goal of [2] to 
relate the relevant generating functions to the Hasse-Weil zeta function of the vacuum 
moduli space. 

One salient and, as it transpires, important feature of these supersymmetric gauge 
theories is that they have an underlying structure of a finite, labeled, directed graph 



^ Quoting P. Candelas, FRS, from his recent London Triangle talk, this mysterious connection, 
which has already produced some beautiful results, beckons for deeper exploration. 
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with a path-algebra encoding the interaction terms in the Lagrangian. Thus was ger- 



minated the industry of quiver gauge theories 21 . In fact, for a vast class of these 
theories, viz., those whose vacuum affords a toric description |26| - consisting of the 
majority of the theories in the AdS/CFT correspondence, as the Calabi-Yau cones 
involved are largely the only ones whose metrics are known explicitly - a simpler still 
graphical description exists. For the myriad of these theories, infinite families thereof 
having already been identified, a bi-partite graph (dimer model) on a torus, or, equiv- 
alently called a periodic brane tiling of the plane, suffices completely to encode the 



physics 27,28 



The study of these graphs has been a vast subject over the past decade and the 
algebraic and differential geometry of the associated Calabi-Yau space, as well as the 
combinatorics emergent from enumerations of inequivalent theories and counting of 
BPS spectra, have been exploited to remarkable and beautiful consequences. It is not 
the intention of the present writing to address these issues. It is, however, our purpose 
to investigate, in the spirit of the number theoretical properties which underly the 
construction. Indeed, in loc. cit., given that the central object is an algebraic variety, 
it was natural to study these affine varieties over finite fields. Here, as another central 
object is a finite graph, it would be most conducive to examine these theories under 
the number theoretical guise of these graphs. 

Luckily, the concept of a zeta function for a finite graph, in the fine tradition of 
Weil et ai, has been long extant: this is the Ihara-Serre zeta function. The subject is 
vast and has been generalized over the decades to encompass all partially directed finite 
graphs. The analogues of primes, functional equations, eigen-spectrum and Riemann 
Hypothesis have all been established or conjectured. 

Thus a natural course of action recommends itself and we shall allow ourselves to be 
freely lead by this path. In ^ we familiarize ourselves with the rudiments of the Ihara 
zeta function and set the nomenclature which will facilitate explicit evaluation. Next, 
in ^ we apply the techniques to the non-chiral quiver theories whose vacuum moduli 
spaces are affine Calabi-Yau varieties. We then proceed to address the bi-partite graphs 
which represent gauge theories which inherently have a toric geometry in ^ Finally, 
we conclude with prospects in ^ 
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2 The Graph Zeta Function 



In his classification of tlie discrete subgroups of PGL{2; k), the projective general 
linear group in dimension 2 over a p-adic number field k, Ihara |3| constructed a zeta 
function of Selberg type; i.e., a zeta function whose analytical structure is determined 
by the spectral data of Laplace operators on a surface. Serre pointed out some years 
thereafter that Ihara's construction had an underlying graphic interpretation ji] and 
thus was born the concept of the graph zeta function. 



In particular, the Ihara-Serre graph is an undirected one and over the years [5j-[9] 
there have been various generalizations and extensions, culminating in the recent paper 



10 which presents a succinct explicit formula for the zeta-function of any finite graph. 



The reader is referred to a thorough introduction to the subject in |11| as well as some 



charming accounts in various lecture notes 12 ; from these we first take some requisite 



rudiments and develop an appreciation before moving onto the aforementioned general 
expressions recently obtained. 

The initial definition of the Ihara zeta function was for graphs G which are finite 
{i.e., having a finite number of edges and nodes, say m and n, respectively), undi- 
rected {i.e., the edges have no direction), connected {i.e., each node can be reached by 
traversing some edge) and usually tail-less by containing no degree-1 vertices {i.e., no 
node which is connected by a single edge only; recall that the degree (or sometimes 
called valency) of a node is the number of vertices adjacent thereto by the connection 
of an edge). These restrictions will be relaxed as we proceed. The graph was, however, 
allowed to have multiple edges and loops. 

The n X n adjacency matrix is such that the {i,j)-th entry is the number of 
edges from node i to j, taken with the convention that the i-th diagonal entry is twice 
the number of self- adjoining loops to the i-th node. We can further assign a standard 
direction in G by making each of the m undirected edges a bi-directional arrow and 
introduce a label ei, . . . , for one set of the arrows and Cm+i = e^^, . . . , = 
the contrary set in the oppose direction. A primitive path P in G is then a closed 
path which never back-tracks or has tails. This is to say that in the path along directed 



edges Cj^ . . . e^^, beginning and terminating on the same node, no cj is ever equal to e^-^^ 
nor can the first arrow be the back-track of the last: e^^ 7^ Cj^. Furthermore, we define 



two closed paths to be equivalent if they are identical but for the choice of the initial 
(and thus terminal) node; the length i{P) of the path P is, naturally, the number of 
edges traversed. At last, we can define the most useful concept for a zeta-function: a 
prime in G is the equivalence class of prime paths. 

It is perhaps helpful to illustrate these concepts diagrammatically and we shall take 
a figure from [IT] (which, for amusement, should be reminiscent of a well-known quiver 
for Hirzebruch, though here undirected) and reproduce in part (a) of Figure flj The 




Figure 1: Examples of finite graphs: (a) a non-chiral quiver whose every edge has been 
replaced by a bi-directional arrow; (b) the Cayley graph for Z3. In our nomenclature, the 
former is a fully directed graph, and the latter, fully undirected. A partially directed graph 
would consist of a mixture of edges and arrows. 

paths [626365], [61626364], and [616263646161064], of lengths respectively 3, 4, and 7, are 
clearly prime, but so is the infinite family [(61626364)^6161064] for any n G Z>i. It must 
be warned that in this path algebra, we do not have unique factorization into primes. 

Thus prepared, the Ihara or Graph Zeta function for a tail-less finite undirected 
graph G is defined to be the formal (possibly) infinite product over primes and in a 
complex variable z: 

Cc{z):= n (2.1) 
[p] prime 

where, we recall, ^([-P]) is the length of the path. Moreover, in the product over 

primes, we distinguish between [P] and [P~^] since orientation reversal is not in the 

same equivalence class, which is reserved, as stated above, only for cyclic relabeling of 

the edges traversed. The similarity to and inspiration drawn from the famous Euler 

product Yl (1 ~ P^^) ^ for the Riemann zeta function C(s) is obvious. The depth 
p prime 
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of this analogy, of course, transcends a mere formal resemblance and the number- 
theoretical groundings have supplanted the motivations for many investigations. Some 
salient results have included the following. 

First, let us give some definitions. We call a (g + l)-regular graph G (z.e., each and 
every node has degree g + 1, and hence g + 1 adjacent neighbours) Ramanujan if the 
maximum of the absolute value of eigenvalues of the adjacency matrix A, excluding q+1 
itself, is bounded by 2y/g. In other words, max{|A| : A G Spec(A), |A| 7^ g + 1} < 
Moreover, the Ihara zeta function of a (g + l)-regular (connected, tail-less) graph is said 
to satisfy the Riemann Hypothesis if for < Re(s) < 1, the zeros of Cg(q'~'*)~^ lie 
on the Re(s) = | line. Armed with such suggestive terminology we have the following 
known results enjoyed by the zeta function for a finite graph of n nodes and m edges 
without tails: 



• Analogue of Weil Conjectures: Cg{.z) is rational and is, in fact, the inverse of a 
polynomial. 

• Analogue of the Prime Number Theorem: defining 7r(m) as the number of primes 
in the graph G of length m, as the greatest common divisor of the lengths of 
primes, Rx the radius of convergence of Cg(^) {i- e-, the position of the nearest pole 
to the origin), then 7r(m) = if Ag does not divide m, otherwise, 7r(m) ~ 

as m tends to infinity. 

• There is a functional equation for the graph xi-function for (g + l)-regular G, 
defined as ici^z) := (1 + ^)--"(l - z^il - qzrCc^z), viz., iciz) = ^dj-J- 

• For (g + l)-regular G, Cg{z) satisfies the Riemann Hypothesis iff G is Ramanujan. 



These above elegant properties are indeed inviting. What is perhaps intimidating 
is the seeming abstruseness in the definition of the Ihara zeta function and a desire for 
a computationally feasible expression has emerged since the original works. Over the 
years, the various authors in cit. ibid, have been writing explicit formulae for various 
generalizations of the class of graphs encountered above. We resort to a recent work 



of 10 which gives a closed form, or a working definition if you will, of the Ihara zeta 
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function for an arbitrary finite graph G: 

(1 _ ^2yTr{Q-I)/2 

^""^^^ ^ det{I-Az + Qz^ + Pz^) ' ^^-^^ 

this shall be the expression which we shall henceforth employ. A few clarifications on 
the notations above: 



• The graph G, with n nodes, is allowed to have both edges (undirected) and arrows 
(directed); such a graph is called partially directed. All subsequent matrices are 
n X n. 

• The matrix A is the undirected adjacency matrix, whose (i,j)-th entry is the 
number of length-one walks, along either edges or arrows, from node i to node 
j. We emphasize that the convention is to have the diagonal entries of A being 
twice the number of self-adjoining loops. 

• The matrix P is the directed adjacency matrix]^ whose {i, j)-th entry is the number 
of arrows from node i to node j. Thus, when G is fully undirected (having no 
arrows) P is the zero matrix and when it is fully directed (having no edges) 
P = A. 

• The matrix Q is the matrix of undirected degrees, defined as follows. For a 
node i, the undirected degree is the number of walks of length 1 starting at i 
using only (undirected) edges. That is to say, it is the number of adjacent nodes 
(including itself) by edges only and not by arrows. Q is then the diagonal matrix 
of the undirected degrees of the vertices subtracted, by convention, by the identity 
matrix. For any graph, one can readily see that the exponent of the numerator 
— Tr((5 — /)/2 is equal to the number of nodes minus the number of edges. Note 
that for any self- adjoining loops, the undirected degree is counted as 2, as is with 
the matrix A, because the beginning and the end of the edge is ambiguous. 



For illustrative purposes, let us give a hypothetical example with the appropriate 

i This matrix distinguishes, for example, the situation where a node has (1) an edge loop or (2)two 
arrow loops, both of which have the same A (note that they are both length 1 because the two arrows 
are in the same equivalence class). In our convention, however, all self-adjoint loops and bi-directional 
arrows will be denoted by edges. 
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matrices indicated 







A 



We see that A gives the full adjacency information of arrows and edges while P encodes 

only the direct arrows. Moreover, — Tr(Q — /)/2 = — ^Tr ^ ^ = as indeed 

there are two nodes with 2 edges (the other being an arrow) so that their difference is 
0. 

Once we have the form of the Ihara zeta function explicitly, we can use the primary 
definition of the graph zeta function as a generating function for enumerating the 
primes, i.e., minimal non-back-tracking (Eulerian) cycles. This is done by taking the 
logarithmic derivative and series expanding 2^ log Cg(-2^)- the n-th coefficient will give 
the number of length-n prime paths. Thus concludes our introduction to the graph 
zeta function and in the next section we will move onto an intriguing category of graphs 
to which we will happily apply our new acquired technology. 



3 Quiver Gauge Theories 

Quiver gauge theories have over the past few decades attracted substantial interest 
both in physics and in mathematics, furnishing, in the former, large classes of gauge 
theories with product gauge groups, of potential phenomenological interest, and, in the 
latter, representations of finite graphs, especially in the context of algebraic geometry. 
The constant dialogue between the two has been most fruitful and we refer the reader. 



for example, to 13 , 14 for an elementary introduction as well as a recent thesis for a 



nice comprehensive discussion 15 



To set nomenclature, our quivers will be finite, labeled, directed graphs, allowing 
for loops, bi-directional arrows and self- adjoining vertices. The corresponding (super- 
symmetric) gauge theory will be taken to have gauge group Y\ SU (Ni) where iVj is 
the integer labeling of the i-th node. The matter content consists of bi-fundamentals 
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so that each arrow from node i to j is a field transforming under the {Ni,Nj) of the 
factor SU{Ni) x SU{Nj) in the gauge group. To fully explore the Ihara zeta function, 
we will adopt the convention that bi-directional arrows as well as arrows which join 
a node to itself will be replaced by an edge so that our quivers could be partially 
directed. A gauge theory which consists only of edges is called a non-chiral theory. In 
representation theory, the vector space C^* will be associated to the i-th node and the 
arrows are maps in hom(C^% C'^^ ). Furthermore, the arrows could obey some (non- 
commutative) relations, whereby generating a path algebra; in physics such constraints 
come from a superpotential which governs the interactions. Let us not, for the present 
purposes of the graph zeta function, dwell upon these; we will return to this point in 



3.1 ADE Quivers and J\f = 2 Gauge Theories 

A natural starting point of our investigations is the ADE quivers which are the Dynkin 
diagrams of the simply laced affine Lie groups. We recall these in Figure [2j As gauge 



1 




Figure 2: The ADE Quivers, i.e., the simply-laced Dynkin diagrams of affine Lie alge- 
bras. The nodes are labeled by dual Coxeter numbers, or dimensions of the irreducible 
representations of the corresponding discrete finite subgroup of SU{2). 

theories, these have gauge groups corresponding to 4-dimensional theories with M = 2 
supersymmetry and the dual Coxeter labels give the SU{Ni) factors. It turns out that 
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there is enough supersymmetry here to complete specify the interactions so we need 
not worry about the superpotentials in any case. 



The famous correspondence of J. McKay 16 states that if one took the fundamental 



two-dimensional defining representation 2 of the finite discrete subgroups F of SU (2) 
- whose classification dates back at least to Klein and which we present below in ( |3.3 ) 
to remind the reader - and formed the Clebsch-Gordan (multiplicity) coefficients aij 

with all the irreducible representations r of F, as in 2 ® rj = ^aijrj, then aij is the 

j 

(undirected) adjacency matrix of the graphs in Figure [2] 



There are beautiful algebro-geometric interpretations of this correspondence wherein 
one looks at the quotient singularity C^/F, all of which can be written as a hypersur- 



face affine variety in C^, whose forms we also present in (3.3). Then, the intersection 
matrix of the exceptional P^-divisors in the Gorenstein resolution of the singularity is 
exactly the aforementioned aij. 



Group 


Name 


Order 


An ~ Z„+i 


Cyclic 


n + 1 


Dn 


Binary Dihedral 


2n 




Binary Tetrahedral 


24 


Er 


Binary Octahedral (Cubic) 


48 


Eg 


Binary Icosahedral (Dodecadedral) 


120 



Singularity 


Hypersurface 


• 


xy + z^~^ = 


Dn-. 


+ y^z + ^""^ = 


Eq : 


x'^ + + = 


Er: 


x"^ + y^ + yz^ = 


Eg : 


x'^ + y^ + z^ = 



(3.3) 

Because C^/(F C SU(2)) are locally Ricci fiat and resolve to Calabi-Yau two-folds 
{i.e., K3-surfaces), the geometrical interpretation lends itself well to string theory. 
Indeed, the world volume theory of a D3-brane transverse to the local Calabi-Yau 
threefold C x ADE is precisely the Af = 2 quiver gauge theory described abov^ 



3.1.1 The Cyclic Series 

Let us compute the Ihara zeta function for the ADE graphs, which are all undirected 
and thus are composed only of edges. As a warm-up, consider the A2 = Z3 case, the 



^Strictly, because of the extra factor of C, to each node in the quiver, there is a doublet of self- 
adjoining loops. We shall consider this situation shortly. 
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quiver is cycle with 3 nodes each of which is connected to the other two by an edge. 
This is, in fact, precisely the Cayley graph for Z3, as drawn in part (b) Therefore, 
A=(l i Yp = 03x3andQ = J3x3,givingus,by(|2:2l),a.W = (l-22;3 + z6)-i, 







which is indeed the reciprocal of a polynomial, as required by the graph analogue of 
the Weil Conjectures. 

We can proceed in general and see that Q is always the identity and P, the zero 
matrix and A, the circulant matrix with near-diagonal I's, whence: 

1-2^ + ^2 ^ = 

CaAz)-'=1 1 n=l (3.4) 

1 - 2^"+! + 22^^+2 n > 2 ; 

the degenerate case of n = is the trivial 1-noded quiver with a single self-adjoining 
loop, for which we need A = Q = {{2}} ,P = {{0}}, giving us (ao{z)^^ = (1 ~" ^)^- 
This geometrically is the flat space C^. For n = 1, we have the Cayley graph for Z2 
which, of course, has no loops, and hence the zeta function is just the identity. 

Indeed, these are all 2-regular graphs, with q = 1, m = n. Recalling our discussions 
above, we can define the graph xi-function as ^a„(-2) = (1 ^ ^)^^"^^''Ca„- For n > 
2, ^A„{z) is simply 1 and trivially satisfies the functional equation and indeed the 
eigenvalues of A, being those of the sum of two roots of unity, are bounded by 2y/q = 2, 
and hence satisfy the Ramanujan property. For n = 0, ^c^{z) is also equal to unity, as 
is the quiver Ramanujan. The situation of n = 1 is degenerate since here q = and 
there is no functional equation to speak of, as too is the Ramanujan property violated. 

Parenthetically, had we adopted the convention of making each (undirected) edge 
a bi-directional arrow, then we would have P = A and Q = —Inxn- This gives us 
(^^_^ = (1 - z2)"(det(l - Az){l - z2))-i = (det(l - Az))~\ which is, in fact, the 
general result for all fully directed graphs |7j. 



Self- Adjoining Loops: we mentioned that one should really add two self-adjoining 
loop to each node in the case of 4-dimensional M = 2 gauge theory, in which case for 
undirected edges, we would simply modify A — )■ A + 2 x 21. The first factor of two is by 
convention and the second is natural both physically and mathematically, furnishing. 
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in the former, the trivial matter content arising from the 1®^ when decomposing the 
fermionic 4 — 1®^ 2 (or bosonic 6 — )■ 1®^ © 2®^ ) from the parent U = A theory 



upon the orbifolding [13 17 , and in the latter, the —2 self-intersection of the exceptional 
divisors in the blowup. From a representation theoretical view-point, the new A is now 
21 subtracted by the Cartan matrix of the corresponding affine Lie algebra. 



Such a modification to the diagonal of A brings about a much more complicated 
form for the Ihara zeta function. A general form is difficult to write down, but we here 
present some first values: 



n 













{z^ 


- 1)^(5^2 -6z + 


1) 




1 


{z - 


- l)\z + lY{hz^ 


-5^ + 1) (5^2 _ 3;2 + 1) 




2 


{z - 


- iy{z + lf{hz- 


1) (5^2 _ 3^ + 




3 


{z - 


- lf{z + lf{hz- 


I) ibz"^ - Az + if {^z"^ - 


2^ + 1) 


4 


{z - 


- l)^\z + lf\hz 


- 1) (25z^ - 35^3 + 21;z2 


-7z + if 


5 


{z - 


-iY\z + iY\hz 


- 1) (5^2 _ 5^ + 1)2 ^^^2 


-3^ + 1)^(5^2-2^ + 1) 



(3.5) 



We note that the addition of the adjoints makes these quivers 4-regular graphs, so that 
g = 3. However, the Ramanujan condition is violated in general because we can see 
explicitly the existence of eigenvalues between the trivial eigenvalue of 4 and the bound 
2v^, for sufficiently large n. 



3.1.2 The Dihedral Series 



Let us now depart the comfortable realm of regular graphs and proceed on to the D and 
E series. With zeros on the diagonals of the A-matrix, we readily find the trivial result 
that for all ra, Cr>„(^) = 1- This is immediate from the fact that there are no cycles in 
the quiver. After all, it is important to have cycles in order to have gauge invariant 

25]). Therefore, we shall include the double self- adjoining loops 
as discussed above by letting the diagonal entries of A be 4, generating much more 



operators (cf. [18 23 
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involved expressions, some incipients of which are: 

n CdM-^ 

1 -{z - + 1)9(2^ - If (28^3 - 16^2 + _ 1) 

2 -{z - + 1)11(2^ - 1)4 (24^3 - 20^2 + 8z - 1) (24^^ - 36^^ + 20;^^ _ 7^ + i) 

3 -{z - lf^{z + lf^{2z - If {2Az^ - 40^3 + 24^^ - 82 + 1) (120^^- 
176^4 + 120^3 - 48;z2 + _ 1) 

-(z - lf^{z + l)^5(2z - If (120^5 - 200^4 + 140^3 - 56^^ ^ - 1) (120;z6 

-272^5 + 260^4 - 148^3 + 52^^ _ Hz + 1) 
-{z - lf^{z + lY\2z - If (120^6 - 296^5 + 300^^ - 172^^ + eOz^ - I2z + 1) (600^^ 
-1360^6 + 1396^5 - 880^4 + 360^^ - 96^^ ^ _ 

(3.6) 

What can one say about the spectral as well as analytic properties? There are some 
partial results for irregular graphs into whose territory we have now entered; one useful 
theorem, due to Kotani-Sunada (cf. Chapter 8 of |11|), is that the radius of convergence 
Rc of Cg(-2) (which we recall to be the position of the pole nearest to the origin) lies 
between the reciprocals of the minimal degree p + 1 and maximal degree g + 1 in G: 

g-' <Rg< p-' . (3.7) 



We see that the radius is rather small indeed, never exceeding unity. In fact, their 
theorem further guarantees that every pole of Cg{z) has absolute value between Rq 
and 1, inclusive. 

On these grounds, the graph theory Riemann Hypothesis is that Cci^) is pole 
free for Rg < < \/Rg, or, in its weak version, for Rq < \z\ < 1/^/^. The motivation 
of this hypothesis comes from the consideration of the substitution ({z := {Rq)'^)- 
Then, we wish for a "pole-free" region inside the "critical strip" < Re(s) < 1, or 
at least, in short of general statements of a functional equation for irregular graphs, 
inside the open half-strip | < Re(s) < 1. It turns out that there is a paucity of graph 
which satisfy this hypothesis, whence the weakened version stated above. In the case 
of regular graphs, where Rq = the strong and weak versions are the same. 

Indeed, our A-series, without the self-adjoints, are regular and Ramanujan (except 
for the degenerate case of n = 1) and do indeed obey the strong Riemann Hypothesis 
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(with Rg as well as all the pole being trivially at 1) while with the addition of the 
double self-adjoints, for n > 5 both versions of the said hypothesis are violated while 
below 5, the strong version is satisfied except for n = 2, which only obeys the weak. 
For the D-series above, both are violated for all n. 

3.1.3 The Exceptional Series 

Finally, for the E-series, we also find, in the absence of self-adjoint loops, that the 
zeta-function is unity because of the lack of cycles. In providing the usual entry of 4 
to the matrix A, more interesting forms for C,g{z)~^ emerge: 

Eq : -{z - lf^{z + If^ {2Qz^ - mz^ + 2Az'^ - Sz + if {l2Qz^ - llQz* + l2Qz^ - A8z^ + llz - 1) 

£^7 : -{z - ly^iz + 1)15 (100^^ - 260^5 + 280^^ - 168z^ + 60z^ - 12z + 1) (2400^^ - 7840^^ 
+11620^^ - 10464^6 ^ 6356^5 - 2704^^ + 804^^ _ iqq^2 ^ _ 

Es : -{z - iy^{z + 1)1^ (1200000^1^ - 8000000^^^ + 25186000^^^ - 49984400^^^+ 

70319040^13 _ 74636224^12 ^ 62013952^ii - 41253888^^° + 22261008^^- 
9801936^^ + 3521088^^ - 1025344^'' + 238784^^ - 43456^^ 
+5952^=^ - 576^2 + 35z - 1) 

(3.8) 

These turn out to not to be special enough to satisfy either the strong or weak versions 
of the graph Riemann Hypothesis; nevertheless we will investigate shortly some spectral 
properties of these and the our previously encountered examples. 

3.1.4 Spectra and Roots 

The concept of the Ramanujan graph was introduced to study the spectral properties, 
i.e., the gap of eigenvalues of the adjacency matrices of regular graphs. The fact, 
as discussed above, that the Ramunujan property and Riemann property of graphs 
are equivalent is indeed elegant. Much work has been focused on generalizing the 
various bounds for the maximal eigenvalues as well as the relation to to the pole-free 
regions of the Ihara zeta function. Along a parallel vein has been some recent interest. 
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continuing again on a rather old path, in the investigation of the roots of constrained 
polynomials (q. v. [19] and especially |20| in the context of gauge theories and Calabi- 



Yau geometries, so too is the reader referred to 31 for interpreting the zeros in the 



critical strip of the classical Riemann zeta function as an eigenvalue distribution). 

It is therefore expedient to conduct some numerical experiments, similar to chapter 
8 of [IT] , of the collective pole positions of the zeta function versus the spectral data of 
the adjacency matrices. In Figure [3| we compare and contrast these roots. In parts (a) 
and (d), we present, respectively, the positions of the poles of the Ihara zeta function 
and the spectrum of the undirected adjacency matrix A (with, as always henceforth 
for the ADE cases, the self- adjoining loops), for the first 100 n in the A^-series. The 
plot (d) is rather simple because here A is symmetric, whereby providing only real 
eigenvalues. In (a), there are two distinctive pole-free regions: one on the real axis 



between —1 and 1/5 (the points ±1, coming from the numerator contribution in (2.2) 
are automatically poles) and the other, the missing conjugate pair of arcs from the semi- 
circle, which does not seem spurious since increasing n does not appear to begin closing 
this gap. In parts (b) and (e), we do the same for the D-series, up to n = 100. Again, 
the symmetry of the adjacency matrix forces real eigenvalues and we densely cover the 
region [0,4]. In (b) a hirsute decoration is added to the semi-circular distribution of 
poles but again the behaviour on the real line is as in the case of (a). Finally, in (c) 
and (f), we plot the same for the exceptional series of -Eej.s- 

As a parting indulgence before we move to other classes of gauge theories, let 
us amuse ourselves with a comparative study. Seeing the forms of the hypersurfaces 



in (3.3), one could not help but be reminiscent of the local Hasse-Weil zeta function, 
which, we recall, is the exponentiated generating function for the number of Fpr-rational 
points of a variety over the finite number field of characteristic q = for some prime p. 
Recently, these functions have been computed for a host of gauge theoretic geometries 
in (2l. 



Of course, because we are dealing with afiine singular varieties, rather than smooth 
projective ones, immediate application of the powers of the Weil Conjectures is not at 
hand. However, for the A-series, having the form xy = z"', the answer is simple. For 
the degenerate case of n = 0, or C^, it is a standard fact that the Hasse-Weil local zeta 
function is Cp{t; C^) = (1 — p'^t)~^. In general, for arbitrary n, we have that at z ^ 0, 
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Figure 3: The distribution of the poles of the lhara zeta function for the A-D-E quivers 
in (a), (b), (c) respectively, with the An and Dn series up to n = 100, as well as the 
eigenvalues of their respective adjacency matrices in (d), (e) and (f). 



for every non-zero x there is a unique y, giving us (g — 1)^ points and for z = 0, the 
equation becomes = xy so that there are 2q — 1 points. In total, there are points 
over ¥q so that the local zeta function remains as (1 — amusingly reminiscent 

of the lhara zeta function at n = 0. Since each of our graphs has an underlying affine 
geometry, it is perhaps fruitful to compare the lhara and Hasse-Weil zeta functions 
systematically, a task we will leave for the future. 



3.2 A Zoo of Chiral Quivers 



The ADE theories which we studied in some detail above is only a very small sample of a 
"zoo" of quiver gauge theories arising from geometry, string theory and particle physics. 
Indeed, due to the fact that they are non-chiral (in that each arrow has a partner in 
the opposite direction between the same two nodes), interests in them have been more 
shown in their beautiful mathematical structure rather than their phenomenological 
importance. There is a vast number of chiral theories such as (supersymmetric) Yang- 
Mills theories or Chern-Simons theories in various dimensions which have product 
gauge groups and bi-fundamental matter, and thereby affording quiver description. 
Many of these have interesting vacuum moduli space prescribed by non-trivial algebraic 
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varieties, so too do many emerge as world- volume theories of branes in string theory 
(cf. [21]). Systematic and step-wise enumeration of these have only recently been 
undertaken (22H24 . 



It is therefore only natural that, familiarized by the ADE examples, we move to 
study the plethora of these chiral quiver theories. The astute reader would point out 
that this is perhaps less meaningful than what was addressed previously, since in the 
ADE case the graph alone determined the physics since the interactions, because of 
the large amount of supersymmetry, can be fixed by the matter content, whereas in 
general we need quite more additional information. In Af = 1 chiral gauge theories, for 
example, we will need to insert a superpotential which in turn imposes an algebra on 
the arrows, or in the case of Chern-Simons quivers, we need levels as well as couplings. 
We will return to this issue in the ensuing section which consists of an infinite family of 
theories completely determined by graphical data; for now, as much as to indulge our 
curiosity as for the sake completeness, let us move onto some representative examples 
of chiral A/" = 1 quivers in 4-dimensions, which generically have Calabi-Yau threefold 
moduli spaces and which have potential phenomenological significance. 

Therefore this brief section will only concern with the quiver part of the physics. Let 
us begin, with a graph quite similar to part (b) of Figure [l| but with 3 uni-directional 
arrows from one node to the next, cyclically; this is the graph for C'^/Zs (otherwise 

known as del Pezzo 0) and can be thought of as the McKay quiver for Z3 in dimension 

/ 3 \ 

three. Exphcitly, it has P = A = I 3 J and Q = Diag({ — 1, —1, —1}), giving us 



C{,z)~^ = 1 — 272;'^. The poles are situated at one-third the three cubic roots of unity, 
therefore the radius of converge is Rg = 1/3 with the entire complex plane being 
pole-free except on the circle of this radius, so indeed it satisfies the strong Riemann 
Hypothesis. 

Other famous 4-dimensional theories include the canonical example of D3-branes 
on C'^, which is a so-called "clover quiver", the precursor to AdS/CFT, or the conifold 
theory which is a SU{N) x SU{N) theory whose classical mesonic moduli space is the 
affine Calabi-Yau singularity uv — zw = 0. In our convention, because we replace all 
bi-directional arrows as well as self-adjoints with a edge, the clover has a single node 



with A = {{6}}, P = {{0}} and Q = {{5}}, so that CgW"^ = (1 - z^i^ - 62 + 5z^), 



the same as the Aq result in the ADE case. The conifold has A = I I I ), P = and 
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Q = Diag({l, 1}), giving us C(^)-^ = (1 - z^. 



As we pointed out for amusement that the first example we ever gave, in part 
(a) of Figure [T| was quite akin to something well-known: the so-called Hirzebruch 
quivers. These are gauge theories arising from D3-branes transverse to an afiine 
Calabi-Yau threefold which is a complex cone over the first Hirzebruch surface, or, 
Fq = X P-*^. We have used the plural because, as it turns out, there are two well- 
known Seiberg dual phases (cf. p33 of |14j), both of which are fully directed, with, 

' 2 ~ 

respectively, P/ = A/ = I "002 1 , = Diag(-1, -1, -1, -1) and P// = 



0002 




0202 
0020 

4 
2 



.Qn = Diag(-1, -1,-1,-1). Whence, we have Ci{z) = (1 - 16^^)-^ 

and Cn{z) = (1 — 32^;^)^^ (indeed both conform to the theorem of Mizuno-Sato [t] 
that for fully directed graphs, ({z) = det(J — zA)~^) and both have poles only on some 
circle and thus obey the strong Riemann Hypothesis. 

Now, let us investigate some partially directed graphs in our present sense. From the 



taxonomy of 24 , we take a few examples of quivers which arise in the study of world- 
volume Yang-Mills theory of D3-branes, or the Chern-Simons theory of M2-branes. 
We have chosen an illustrative admixture of partially directed graphs, having various 
nodes, or combinations of self-adjoining loops and bi-directional arrows represented 
by edges. The quiver in (a) has been affectionately dubbed the "phallus" and has 
recently been of interest in the study of M2-branes. We see that the relevant matrices 
are v4 = ° P = and Q = Diag({0,4}), giving us ({z)-^ = (1 - z^){l - 

4z + 3z^). The theory in (b) is the graph for the so-called "suspended pinched point" 
(SPP), an affine Calabi-Yau threefold variety described by xy = zw"^; it is a theory 

which dates back to the early days of D-brane technology. Here we have, as a slight 

/ 2 1 1 \ 

modification of the Cayley for Z3, ^ = I 1 » 1 , P = and Q = Diag({3, 1, 1}). 

V 1 1 y 

This gives us ((z)^^ = —{z — 1)^ (3^;^ -|- Az^ + Az'^ + Iz^ — 1). Finally, the quiver in 
/ 1 \ / 1 \ 

(c)hasA= h; : 2 ; \\ andg = Diag({-l,-l,0,0}). This 





gives us = (1 — -2^ — 2^;^ + -z^) ^ ■ All the three Ihara zeta functions in this set of 
examples are seen to satisfy the strong Riemann Hypothesis. 



19 




Figure 4: Some samples of quivers which have arisen in the study of D-brane gauge theories 
and M-brane Chern-Simons theories. These are all partially directed graphs, having either 
self-adjoining loops or bi-directional arrows represented by edges. 



4 Dimer Models 



We raised the objection in the previous section, after leaving the comfortable terrain of 
ADE-quivers, that an immediate trouble with the zoo of quivers arising from general 
Yang-Mills or Chern-Simons theories is the insufficiency of the graph to completely 
specify the physics and the extra information requisite seem to lie beyond the scope 
of the zeta function. Nevertheless, there is a vast number of tremendously interesting 
theories which have of late generated much beautiful mathematics and physics alike 
(cf, reviews in [27|[28]). These are completely specified by a finite graph, viz., a bi- 
partite graph (a dimer) drawn on an torus, or, equivalently, a periodic brane tiling of 



the plane, and has and underlying affine toric variety as the vacuum moduli space 26 



We will not delve into the immensity of the subject here but will recall some key 
definitions first. A bi-partite graph is one whose set of nodes can be partitions into two 
disjoint sets B and W such that no node in B is adjacent to any other nodes in B and 
likewise for W. Graphically, we may represent the nodes in B as, say, "black", and 
those in W, as "white", so that our dichromatic graph will have edges only between 
white and black nodes. 

Consider such a graph as doubly periodically extended on a plane {i.e., on a torus), 
with a finite representative in the fundamental region F. Letting each face formed by 
loops in the edges in F represent a gauge group factor, each edge, a matter field, such 
that a black (respectively white) node of degre^ k contributes a degree k monomial 
term with a plus (respectively minus) sign. Then our bi-partite graph represents a 



^Or sometimes called valency. 
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quiver gauge theory, a construction which has come to be dubbed dimer model or 
brane tihng. It is an elegant and important fact that the majority of (superconformal) 
gauge theories known to the AdS/CFT Correspondence fall into this class, being thus 
representable. 

For the sake of computing the zeta function, we will not need anything as refined as 
the Kasteleyn matrix, which is a weighted adjacency matrix who determinant counts 
the number of perfect dimer matchings as well as the being intimately related to the 
Newton polynomial of the underlying toric variety which describes the vacuum moduli 



space of the gauge theory 27-29 . The matrices in (2.2) shall suffice for our present 
purposes. 



4.1 Tilings, Bi-Partite Graphs and Ihara Zeta 



To the beautifully compiled Appendix A.l of 22 we shall henceforth refer the reader, 
for the notation, the corresponding gauge theory and, especially, the periodic tiling by 
the bi-partite graph. The programme, undertaken in the aforewritten, of extracting the 
various adjacency and degree matrices and computing the Ihara zeta functions, lends 
itself easily to us. Indeed, because the edges in the dimer have no preferred direction, 
in the following examples all our P matrices will be zero. 

Let us commence with the simplest tiling, corresponding to the "clover" quiver of 
the TV = 4 Yang-Mills which launched AdS/CFT and whose vacuum moduli space is 
the trivial affine Calabi-Yau threefold C^. The dimer is a hexagonal periodic tiling, 
with one pair of black/white nodes, each of undirected degree (valency) 3, given in (1.1 
of A.l in cit. ibid., which, for illustration, we include in Figure 
and Q = Diag({2,2}), giving us C(-2)"^ = (1 -'2^)(1 - 5z + Az"^ 
these matrices is twice that of our familiar Kasteleyn matrices, this is simply because 
we are not addressing issues such as perfect matching and are concerned only with the 
adjacency information; hence a pair of black- white nodes requires dimension 2 matrices 
rather than 1. 

It is easy to see, that in all our cases, A = K ^ j with K the unweighted 

Kasteleyn matrix because we, by our convention, use the absolute adjacency and not 



Hence, A .30 
Note that the size of 



3 
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only with the rows indexing one colour and the columns, the other. Here, is the 
standard Kronecker tensor product of matrices. Similarly, Q = (^l ° j — / , where 
R is the (undirected) degree matrix encoding the valency of the nodes. The reader is 



referred to the interesting paper 30 which has recently cast R in terms of ramification 



indices of Groethendieck's dessin d'enfants in mapping, in the sense of Belyi, the torus 
on which the bi-partite graph is drawn, to a P^. In fact, further simplifications ensue. 



Because of the so-called "toric condition" dubbed in 26 , each matter field appears 



precisely twice in the superpotential, with opposite sign, our bi-partite graph has the 
further property that no white (black) node is adjacent to more than one black (white) 
node. Thus the matrix K is actually diagonal, and is simply R. In summary, for all 
our dimer models, say of n pairs of black-white nodes. 



where R = Diag(rj) has, as its diagonal entries, the degrees of the nodes. 



(4.9) 



c: 




Figure 5: The dimer model for A/" = 4 super- Yang-Mills. The three figures represent, from 
left to right, (a) the quiver diagram for the matter content, the single node corresponding 
to the U{N) gauge group and the three self-adjoint loops, the three adjoint fields; (b) 
the hexagonal periodic bi-partite graph tiling of the plane and (c) the toric diagram of C'^, 
the underlying Calabi-Yau geometry. For completeness, the superpotential is also included 
underneath. The figure is drawn by the marvelous Mathematica package written by J. 
Pasukonis. 



Aided by (4.9), we move on. Next, the conifold theory is a square tiling, with 
one pair of nodes so that A = (^^ ^ j and Q = Diag({3,3}), yielding C(-2)~^ = 
(1 — 2;^)^ {9z^ — lOz'^ + 1). Returning to the subject of {q + l)-regular graphs, it is a 
theorem, following Perron- Frobenius, that — (g-|-l) is in the spectrum of the adjacency 
matrix A iff the graph is bi-partite. Out above two examples are, respectively, bi-partite 
3- and 4-regular graphs, and indeed —3 and —4 are eigenvalues. 
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Moving on to (C^/Z2) x C, presented in (2.1) of Appendix A. 2, Loc. Cit.. The quiver 
for this is actually that of the Ai singularity addressed earlier; however, for the dimer, 

3 

we have two pairs of nodes such that ^ = ( 3 " o o 1 and Q = Diag({2, 2, 2, 2}), 



3 
3 
3 

-1 — n _ ^2\2 /if^^S _ /in-v6 I QQ^4 _ in-v2 



and thus ({z)-^ = (1 - ^ ) (16;^^ - mz^ + 33;^^ - IQz^ + 1). 



The next example is the SPP theory by which we casually passed in ^3.2; now let us 



study it more systematically. The bi-partite graph is given in (2.2) of their Appendix 

(0 3 \ 
3 ° ° g I and Q = Diag({2, 3, 2, 3}), 
4 / 

SO that C(^)~^ = (1 ~ (362;® — 85^;^ + 632;^ — 152;^ + 1). Amidst these above exam- 
ples, the Ihara zeta functions all satisfy the strong Riemann Hypothesis, except SPP, 
which violates even the weak version. 

All periodic tilings by bi-partite graphs up to 8 pairs of black/white nodes, which 
are consistent in the sense of being well-defined gauge theories, have been classified in 



the lovely work of 22 and a computer database thereof, compiled. Our examples above 
comprise but an hors d'oevres to this delicious database which, perhaps surprisingly - 
in contrast to the expectation of a myriad of contingencies - boasts only 41 consistent, 
inequivalent members. This number is small enough to afford us presenting them all 
in Appendix |A| In it, we present the quiver and tiling data and compute the zeta func- 
tion associated to these finite graphs, indicating, in particular, whether strong/weak 
versions of the graph theoretic Riemann Hypothesis are obeyed or violated. For amuse- 
ment, as in §3.2[ since we are naturally confronted, in addition to our bi-partite graphs, 
with a set of (in general) chiral quivers, it is immediate that we should compute the 
zeta function for these graphs as well, a simple exercise we do undertake and present 
together with our aforementioned appendix. 



As was discussed in ^3.1, a comparative study of the pole-free region and the spec- 
trum of the full adjacency matrix is illuminating. We perform this task and give the 
results in Figure [6j In part (a), a compounded plot of the poles of the Ihara zeta 
function (marked as "P") and the spectrum (marked as "E") of the adjacency ma- 
trix (double the size of the unweighted Kasteleyn matrix) is drawn; to the simplistic 
appearance of the plot we will shortly return to remark. In part (b) and (c), we respec- 
tively conglomerate the poles of the zeta functions and the eigenvalues of the adjacency 
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matrices for all the quivers associated to the 41 consistent tilings. 

Jit- . 




-2 

-2.24 

-s,s - 

-J (b) 

Figure 6: (a) The combined plot of the poles of the lhara zeta function (marked "P") and 
the spectrum (marked "E" ) of the adjacency matrix for the bi-partite graphs corresponding 
to the 41 consistent dimer models classified in the literature, (b) The zeta poles and (c) 
the spectrum of the adjacency matrix of the associated (generically chiral) quivers. 




The simplicity of part (a) of the figure is of no surprise to us. From (2.2) and 



(4.9), and recalling that — Tt{Q — I)/2 is equal to the difference between the number 
of nodes and edges, we can readily obtain the analytic expression of the zeta func- 
tion. Letting the list of valencies (ramifications) of the dimer be {ri,...,r„}, the 

n 

difference between nodes and edges is simply 2n — ^r^, which given that each rj 

1=1 

is at least unity, will always be less than 0, whereby bringing the numerator down 
to the denominator and ensuring that the zeta function remains the reciprocal of a 
polynomial, as required. Now, because the denominator has P = 0, it reduces to 



det 



Diag (1 + (ri - 1)^2, . . . , 1 + (r, - l)z^) ® {ri, rs, . . . , r„} ® 



1 

1 



givmg us. 



upon expanding the determinant of the block matrix. 



2n-\- r ■ 

C(;.)-^ = (1 - z') - ' Hil + (r, - l)zY - {r^zf . (4.10) 



=1 



Thus, the roots are at ±1 and ±(1 — r^) ^, all real and between —1 and 1, as indicated 
by the figure. Now, the radius of convergence is the distance of the closest pole to 



the origin, thus, Rg = {rmax — 1) \ the integers p, g in (3.7) are p = Tmin — 1 and 
Q. = fmax — 1, so indeed Rq G [q^^jV^^] and saturating the lower bound, so that the 
strong and weak versions of the Riemann Hypothesis coincide, which here state that 
the region {rmax — 1)^^ < l-^l < ('"max — 1)~^ is pole free in the complex plane. In 
other words, all must be that they must satisfy the inequality: either (r^ — 1)^^ > 
{'fmax — 1)~^ or (rj — 1)^^ < {vmax — 1)^^- The second inequality being vacuous, we 
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are thus lead to the conclusion that a dimer model/brane-tiling satisfies the Riemann 
Hypothesis (as mentioned above, here the strong and weak versions equate) iff 



2ri + 2 < r„ 



(4.11) 



for all valencies (ramifications) less than the maximum. 



It is entertaining to meddle with a little statistics. In 22 , the state of the art 
classification of dimer models was carried, including even the inconsistent ones in a 
stringent sense that they have good IR behaviour. Nevertheless, relaxing this condition, 
we are granted a database of 375 graphs into which one could liberally dip. Of these, 
it is found that 112, or about 40.1% obey the Riemann Hypothesis. We also repeat 
Figure |6] for these members, and not merely the distinguished 41. 
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Figure 7: (a) The combined plot of the poles of the lhara zeta function (marked "P") and 
the spectrum (marked "E" ) of the adjacency matrix for the bi-partite graphs corresponding 
to all the known dimers, up to valency 8 and not necessarily consistent gauge theories, (b) 
The zeta poles and (c) the spectrum of the adjacency matrix of the associated (generically 
chiral) quivers. 



5 Conclusions and Prospectus 

We have taken an adventurous journey accompanied by two seemingly distant acquain- 
tances, one belonging to the world of supersymmetric gauge theories with continuous 
vacua and the other, to zeta functions with analogues of functional equations and zero- 
free regions. The common language is that of finite graphs, either in the form of quiver 
representations or in the guise of periodic bi-partite tilings of the plane. 

Guided by some recent explicit formulae which generalize the classical result of 
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Ihara and which are apphcable to any partially directed graph containing bi-directional 
arrows, cycles or self- adjoining loops, our journey is inevitably driven by Wanderlust 
and experimental urge. Given the extensive catalogue, compiled over the past decade, 
of a myriad of gauge theories with various supersymmetry and underlying algebraic 
geometry, we ventured to compute the Ihara zeta-function for a host of examples. 

We began by studying theories in 3 + 1-dimensions with M = 2 supersymme- 
try, whose moduli spaces are the local K3 surfaces, i.e., quotient singularities of C^. 
The graphs for these are the famous ADE quivers of McKay with two additional self- 
adjoining loops per node. We find the curious fact that only the A-series for n < 5 
satisfy the graph theoretical Riemann Hypothesis (with n = 1 satisfying merely the 
weak and the others, the strong), whilst all the other quivers do not have interesting 
pole-free regions and violate the conditions requisite for the hypothesis. 

We proceeded to the zoo of quivers of less supersymmetry, computed the Ihara 
zeta functions and studied their properties. Then, anticipating the reader's apprehen- 
sion that in these cases the quiver alone does not capture the physical data and that 
formal relations coming from the superpotential need to be imposed on the arrows, 
we focused primarily on the theories whose vacua are affine toric varieties. Not only 
are such theories central to the AdS/CFT correspondence, they are also completely 
depicted by a bi-partite periodic tiling of the plane. We investigated the Ihara zeta 
function for these dimer graphs and showed how to express the Graph Riemann Hy- 
pothesis in terms of inequalities in the valency numbers (order of monomials in the 
superpotential) of the black and white nodes. In due course, we also examined, as 
inspired by Ramanujan graphs, the eigen-spectrum of the full adjacency matrices as 
collected roots of constrained polynomials. 

It is our belief and hope that our excursion had not been a mere self-indulgence but, 
rather, hinted at an amusing diversion on gauge theories, leaving us with the after-taste 
of many a question. What physical quantity does the Ihara zeta function measure 
or represent? How do theories whose underlying graphs satisfy the strong or weak 
versions of the graph Riemann Hypothesis distinguish themselves? What structure to 
the space of conglomerate eigenvalues of the adjacency for the set of consistent quiver 
gauge theories could one uncover? To these and other enquiries let us turn, but on 
another day. 



26 



Acknowledgements 



Ad Catharinae Sanctae Alexandriae et ad Maiorem Dei Gloriam. We would like to 
extend our gratitude to A. Lauder and J. Pasukonis for valuable discussions and further 
to the gracious patronage of the Science and Technology Facilities Council, UK, for 
an Advanced Fellowship, the Chinese Ministry of Education, for a Chang- Jiang Chair 
Professorship at NanKai University, the National Science Foundation, USA, for grant 
CCF-1048082, as well as City University, London and Merton College, Oxford, for their 
enduring support. Above all, to the happy occasion of the author's engagement to the 
lovely Miss E. K. Hunter this humble note is dedicated. 



A Ihara Zeta Functions for Brane Tilings 

In this appendix, we catalogue the Ihara zeta functions for the 41 consistent brane- 



tilings (dimer model) classified in the literature 22 , each of which is a bi-partite 
graph drawn on a torus, and to which is associated a quiver gauge theory of A/" = 1 
supersymmetry and living in 3 + 1-dimensions (generalization to 2 + 1-dimensional 
Chern-Simons theories in also straight-forward). 

The first column is the full adjacency matrix of the quiver diagram for the theory. 
The second column records the list of the valencies or in/out degrees, i.e., the number 
of edges emanating from the black/white nodes. These two matrices completely specify 
the tiling and the reader is referred to [22] for the elegantly drawn graphs. We then 
compute, in the third column, the zeta function for this information (noting that in- 
deed, they are all reciprocals of polynomials) and check, in the fourth column, whether 
they obey the strong ("S"), weak ("W") versions of the graph Riemann Hypothesis, or 
violate it ("N"). For amusement and completeness, we also compute the zeta function 
for the quiver itself (which only partially encodes the physical data) and check their 
relation to the Riemann Hypothesis. 
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